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Traversible wormhole space-times are found as static, 
spherically symmetric solutions to the Einstein equations with 
ingoing and outgoing pure ghost radiation, i.e. pure radia- 
tion with negative energy density. Switching off the radiation 
causes the wormhole to collapse to a Schwarzschild black hole. 
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The idea of space-time wormholes originated with 
Wheeler [jjj, with traversible wormholes becoming popu- 
lar after the article of Morris & Thorne . More recently, 
renewed interest has focussed on outstanding questions 
such as how to construct wormholes || , their dynamical 
behaviour [Q, their stability an d the nature of the 
negative-energy matter needed to support them Q . This 
letter proposes a very simple model of negative-energy ra- 
diation and exhibits resulting wormhole solutions, which 
can be used to study the other questions. 

The simplest radiation model is pure radiation, also 
known as null dust, for which the energy tensor is pu ® 
u, where u is a null vector and p is the energy density 
(which, if positive, can be absorbed in the normalization 
of u). This represents incoherent radiation and occurs 
in the geometric-optics limit of massless particles such as 
photons. The matter model proposed here is simply pure 
radiation with negative energy density. In quantum field 
theory, this is associated with ghost fields, so it may be 
called pure ghost radiation. One might expect it to be a 
similarly valid model of negative-energy radiation. More 
realistic models, like semi-classical quantum fields such as 
those produced by Hawking radiation |^|, are often not 
analytically tractable, so a simple model has theoretical 
merit. 

With pure radiation of the usual positive energy den- 
sity, the static, spherically symmetric Einstein system 
has been studied by Date ||, with solutions found by 
Kramer |ll]] and, more generally, by Gergely Ell]. Re- 
versing appropriate signs leads to the wormhole solutions 
derived below. However, the geometry is quite differ- 
ent, in particular being non-singular, ft is curious that 
benign wormhole solutions naturally appear under such 
situations. 

Liberally adopting the method and notation of 
Gergely, the line element of a static, spherically sym- 
metric space-time may be written as 

ds 2 = r 2 dfl 2 + fir^dr 2 - h{r)dt 2 (1) 

where r is the areal radius, t is the time, dfl 2 refers to 



the unit sphere and the metric functions / and h are 
positive. The active gravitational mass-energy (lj,[ll| of 
such a space-time is 



i(r) = (1 - /(r))r/2. 



(2) 



The energy tensor of pure ghost radiation is —tu®u, 
where u is a null vector and r > is the negative energy 
density. With both ingoing and outgoing radiation, the 
energy tensor is 



T = —T + U + ® U + — T-U^ ® U- 



(3) 



where u± are null vectors. For the above metric, one may 
take 



(4) 



so that the null vectors are relatively normalized. Static 
solutions require the same negative energy density t = t± 
for both ingoing and outgoing radiation. This specialized 
matter model can also be regarded as an anisotropic fluid 
with density — r, radial pressure — r and vanishing tan- 
gential pressure. Here and throughout, units are such 
that the speed of light and Newton's constant are unity. 

The energy-momentum conservation equations V • T = 
reduce to 



(hr 2 







(5) 



where the prime denotes dj dr. Thus the (gravitationally 
red-shifted) radial tension or negative linear mass density 



A = Airr 2 hT 



(6) 



is a positive constant. It is possible to absorb the mag- 
nitude of A in the static Killing vector, but it will be 
retained here. In terms of the function 



= -2X/h 

the Einstein equations G — 8irT reduce to 

rf = 1 - P - f 
rfh> = h(l + (3-f) 
rf/3' = -/3(l + /?-/) 



(7) 



(8) 
(9) 
(10) 



where one equation is redundant, having already solved 
the conservation equation. Transforming to p — lnr and 
using a dot to denote dj dp, the reduced system is 



1 



(11) 

(12) 



(13) 



f=l-P-f 
//? = -/3(l + /?-/). 

The second equation yields 

which can be eliminated from the first equation to yield 

(1 + 0)00 -{[3 + 2)f3 2 + 0(1 - 2/3)0 + 20 3 = 0. (14) 

Dividing by 3 , this becomes an exact differential, inte- 
grating to 



l + (3 



1 







(15) 



where D is the integration constant. This leads to the 
algebraic relation 



(l + /3) s 
2//3 



ln( — 5 

a(3 



(16) 



where a = exp(— (1 + D)/2) is a positive constant. The 
branch of the logarithm corresponds to negative (3, re- 
flecting the negative energy density. Returning to the 
reduced system, one may eliminate p or r to give 



df _ /(/? + /-!) 
d(3 0(0 -/+!)' 



(17) 



Since solutions with negative (3 are being sought, the 
method of Gergely is now modified by defining 



V - 



l + (3 



Then the above equation becomes 



(18) 



(19) 



Integrating the corresponding homogeneous equation 
yields 



p = -2e l (f> 
and the solution is completed by 



^(l) = / e- £ d£ + b= ^—ei£(l) + b 
Jo 2 



(21) 



where b is constant and erf denotes the error function, the 
symmetric integral of the standard normal distribution. 
The solution is mathematically familiar compared with 
the case of positive energy density. Then 



f3=pl-l 



-1 - 2le l 



(22) 



Collecting results, the class of solutions depends on the 
parameters (a, b) and the normalizable A, and may be 
summarized by 



+ 21(f)) (23) 
(24) 

(25) 
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l + 2le l2 (f) 



as substituted into the line element ([!]). By construction, 
this is the unique class of static, spherically symmetric 
solutions to the Einstein equations with pure ghost radi- 
ation. 

It is straightforward to calculate 



2a(f) 
-i 2 



+ 21(f) -2e l (f> 2 )a/2 



dr 

n 

m 

dm / ,2 , 

— = - [l + 2le l 0) acf). 



(26) 
(27) 
(28) 



The line element may then be written explicitly in (t, I) 
coordinates as 



ds 2 = a 2 (er l2 + 21(f)) dn 2 + 2a 2 e~ l2 (e^ 2 + 2l<f>) 



2Xdt 2 



1 + 2le l \ 



2l(f>) dl 2 
(29) 



In the case b = 0, (f)(1) is an odd function and the met- 
ric is even in the spatial coordinate I. This describes a 
symmetric wormhole with spatial topology R x S 2 and 
minimal surfaces at the wormhole throat I — 0, with ra- 
dius r — a. The space-time is not asymptotically flat, but 
otherwise constitutes a Morris- Thorne wormhole. There 
are no singularities, unlike in the corresponding positive- 
energy solutions. The 6^0 cases include asymmetric 
wormholes which are analogous to the asymmetric Ellis 
wormholes E3 for a ghost Klein-Gordon field. 
The solutions may be written in dual-null form 



ds 2 



2 dfl 2 -hdx + da 



(20) 

in terms of null coordinates x ± defined by 



dx~ 



a / 



dt±—=[e 



21(f)) dl 



(30) 



(31) 



Integration by parts yields an analytic solution: 

(l + 2l 2 )(f>). (32) 



x ± =t± 



2^X V 



Thus the metric functions are implicitly known as func- 
tions of x + — x~ . Using the relations 



d±l = ± 



2a(e~ 



21(f)) 



(33) 



2 



the solutions have been checked by substitution in the 
dual- null form of the Einstein equations (l^l^l for the 
general spherically symmetric line element ([50 ) with non- 
zero energy components T±± = —hr/2. 

Returning to more general issues, the simplicity of pure 
ghost radiation means that it is easier to study how the 
wormholes react to changes in the radiation level. This 
has not been analytically tractable previously except in 
a two-dimensional model In particular, the ingoing 
and outgoing radiation simply pass through one another 
without interaction, following null geodesies with propa- 
gation equations d±(hr 2 T±) = 0. For instance, it is not 
difficult to see what happens if the radiation supporting 
the wormhole is turned off, as follows. 

Consider a static wormhole for {x + < 0, x~ < 0}, with 
the ghost radiation then switched off from both sides of 
the wormhole, so that r T = for x > 0, respectively. 
The solutions in the regio ns j x^ > 0,x T < 0} will be 
similar to Vaidya solutions pa] , while the solution in the 
future region {x + > 0, x~ > 0} will be vacuum and there- 
fore, by a version of Birkhoff's theorem valid through 
the horizons of a black hole p2| , Schwarzschild. Fur- 
ther, this region can be seen to be be the interior of a 
Schwarzschild black hole with mass a/2, by differentia- 
bility of the metric at x + = x~ = 0. Specifically, conti- 
nuity of the trapping horizons d±r = |]l2| , |l3i , which 
initially form the throat {x + = x^,x^ < 0} of the 
wormhole, means that they can be joined to the event 
horizons {x^ = 0, x T > 0} of a black hole. This gen- 
erally unexpected behaviour was predicted by a general 
theory of wormholes and black holes p I and recently con- 
firmed in the two-dimensional model 1 3J . In summary, if 
the supporting ghost radiation is switched off, the worm- 
hole collapses to a Schwarzschild black hole. Details of 
such collapses and other dynamical processes involving 
the wormholes are to be presented in forthcoming work. 

Research supported by Korea Research Foundation grant 
KRF-2001-015-DP0095. Thanks to Sung- Won Kim for 
support and Hisa-aki Shinkai for discussions. 



[11] L A Gergely, Phys. Rev. D58, 084030 (1998). 

[12] S A Hayward, Phys. Rev. D53, 1938 (1996). 

[13] S A Hayward, Class. Quantum Grav. 15, 3147 (1998). 

[14] H G Ellis, J. Math. Phys. 14, 395 (1973). 

[15] P C Vaidya, Proc. Indian Acad. Sci. A33, 264 (1951). 



[1] 
[2] 
[3] 

[4] 
[5] 
[6] 
[7] 

[8] 

[9] 
[10] 



J A Wheeler, Ann. Phys. 2, 604 (1957). 
M S Morris & K S Thorne, Am. J. Phys. 56, 395 (1988). 
S A Hayward, S-W Kim & H Lee, Phys. Rev. D65, 
064003 (2002). 

S A Hayward, Int. J. Mod. Phys. D8 . 373 (1999). 



C Armendariz-Picon, ^r-qc/0201027 



K A Bronnikov & S Grinyok, |gr-qc/0201083 . 
N Dadhich, S Kar, S Mukherji & M Visser, Phys. 
D65, 064004 (2002). 

D Hochberg, A Popov & S V Sushkov, Phys. Rev. 
78, 2050 (1997). 

G Date, Gen. Rel. Grav. 29, 953 (1997). 

D Kramer, Class. Quantum Grav. 15, L31 (1998). 



Rev. 



Lett. 



3 



